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We present a combined experimental and theoretical study on strong-field ionization of a three- 
dimensionally oriented asymmetric top molecule, benzonitrile (C7H5N), by circularly polarized, 
nonresonant femtosecond laser pulses. Prior to the interaction with the strong field, the molecules 
are quantum- state selected using a deflector, and 3-dimensionally (3D) aligned and oriented adi- 
abatically using an elliptically polarized laser pulse in combination with a static electric field. A 
characteristic splitting in the molecular frame photoelectron momentum distribution reveals the po- 
sition of the nodal planes of the molecular orbitals from which ionization occurs. The experimental 
results are supported by a theoretical tunneling model that includes and quantifies the splitting in 
the momentum distribution. The focus of the present article is to understand strong-eld ionization 
from 3D-oriented asymmetric top molecules, in particular the suppression of electron emission in 
nodal planes of molecular orbitals. In the preceding article [Dimitrovski et al., Phys. Rev. A 83, 
023405 (2011)] the focus is to understand the strong-eld ionization of one-dimensionally-oriented 
polar molecules, in particular asymmetries in the emission direction of the photoelectrons. 

PACS numbers: 33.80.Rv, 33.80.Eh, 42.50.Hz, 37.20. +j, 37.10.Vz 
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I. INTRODUCTION 



Photoelectron spectroscopy has been and remains a 
ubiquitous and useful method to obtain information 
about the electronic structure of molecules. In the most 
traditional form absorption of a single photon frees a sin- 
gle electron from the molecule and the kinetic energy is 
recorded. The classical use of photoelectron spectroscopy 
is the determination of binding energies of molecular or- 
bitals and possibly of the vibrational substructure [1]. 
Over the past 10-15 years such spectroscopy in combi- 
nation with femtosecond pump-probe methods has been 
demonstrated to be also a powerful technique for obtain- 
ing time-resolved insight into the coupled electronic and 
vibrational dynamics occurring after photoabsorption in 
polyatomic molecules [U [3] . 

In addition to the kinetic energy, the emission direction 
of the electron is a highly useful experimental observable 
provided it is detected for defined vibrational and rota- 
tional states of the cation formed [4] or if it is measured 
with respect to the molecular frame p^ IS] . One experi- 
mental approach to obtaining molecular frame photoelec- 
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tron angular distributions (MFPADs) is to determine the 
molecular frame after the ionization event. This is pos- 
sible if ionization leads to fragmentation of the molecule 
and the electrons can be measured in coincidence with 
fragment ions recoiling in directions that directly reflect 
the molecular orientation at the moment of ionization, a 
condition often referred to as the axial recoil limit. Such 
coincidence methods have been applied successfully, pre- 
dominantly to small molecules where synchrotron radi- 
ation removes an inner shell electron [6]-[9], but also to 
situations where valence electrons are removed by UV 
radiation [10Hl2]. 

The other experimental approach to MFPADs is to 
fix the orientation prior to the ionization event, which 
offers the obvious advantage that no fragmentation of 
the molecule is needed. Using resonant photoexcitation 
mild ID alignment can be created for molecules in ex- 
cited states O [13]. Alternatively, alignment methods, 
based on the application of moderately intense, nonreso- 
nant laser pulses, can confine the principal axes of polar- 
izability of molecules sharply along space-fixed axes [M] 
and, thereby, provide the desired fixed-in-space targets 
for photoionization of molecules possessing an inversion 
center. The first experiments showing MFPADs mea- 
surement from aligned, inversion symmetric (unpolar) 
linear molecules were reported 2008-2009 [llHII]. Tar- 
gets of molecules lacking an inversion center (i.e., polar 
molecules) must, in addition to having their axes con- 



fined, possess a directional order of their dipole moment, 
i.e., the molecules must be oriented as well as aligned 
[14] . One approach to efficient orientation and alignment 
is through the combined action of a laser pulse and a 
weak static electric field [ T8H2T] . 

In the present work we study single ionization of ID 
and 3D aligned and oriented samples of the prototypical 
asymmetric top molecule benzonitrile by intense, circu- 
larly polarized femtosecond laser pulses, extending our 
first results reported recently [22 . For ID aligned and 
oriented molecules pronounced up-down asymmetries of 
the electron emission are observed resulting from a higher 
probability for ionization when the circularly polarized 
field points along the permament dipole moment of ben- 
zonitrile compared to the opposite direction. These re- 
sults are similar to those obtained for the linear OCS 
molecules reported in the accompanying paper [23]. For 
3D aligned and oriented molecules striking suppression of 
electron emission in the polarization plane of the ionizing 
laser is observed. Our theoretical analysis, relying on a 
modified tunnel ionization theory, identifies the suppres- 
sion as resulting from nodal planes in the orbitals from 
which electrons are removed. Note that effects of nodal 
plane structure and orbitals below the highest occupied 
molecular orbital on strong field ionization of in sym- 
metric molecules was discussed theoretically some time 
ago [Ml [25]. 

The paper is organized as follows. In Sec. II, the exper- 
imental setup is discussed. The experimental results are 
presented in Sec. Ill, including a discussion of the method 
of deflection, alignment and orientation of the benzoni- 
trile molecules, and the MFPADs from strong-field single 
ionization of benzonitrile by a circularly polarized fem- 
tosecond laser field. Section IV contains the theory in- 
cluding a discussion of the characteristic off-nodal-plane 
emission pattern. In Sec. V, the theoretical MFPADs 
are compared with the experimental ones. Section VI 
concludes and gives a brief outlook. The molecular prop- 
erties of benzonitrile, used in the theoretical model, are 
summarized in Appendix A. 




FIG. 1: (Color online) Exploded view of the molecular beam 
machine consisting of (from the right) the source chamber, the 
deflector chamber and the target chamber. In this work two 
laser beams are used, the 1064 nm YAG beam (blue cylinder) 
and the 800 nm probe beam (red cylinder). An additional 
800 nm femtosecond beam (green cylinder) can be included to 
induce, for instance, rotational dynamics. The pair of circular 
discs in each laser beam represents a half- wave and a quarter- 
wave plate used to control the polarization state of the laser 
pulses. 



cold molecular beam without cluster formation. We es- 
timate the rotational temperature of the molecular gas 
with two independent techniques. One relies on compar- 
ing experimentally observed degrees of alignment with 
theoretical calculations for a temperature averaged sam- 
ple of molecules [28 , whereas in the other method verti- 
cal molecular beam profiles of electrostatically deflected 
molecules are measured and compared with simulations 
[29^. We find that the typical rotational temperature is 
~1 K. The experiments are performed at 20 Hz limited 
by the repetition rate of the alignment laser. The molec- 
ular beam is skimmed twice, 14 and 38 cm downstream 
from the nozzle of the Even-Lavie valve, respectively, to 
collimate the expansion. The first skimmer, positioned 
at the exit of the source chamber, has a 3-mm diameter 
(Beam Dynamics model 50.8). In the deflector chamber 
the molecular beam passes through the second skimmer 
(1-mm diameter Beam Dynamics model 2) just before 
entering the 15-cm-long electrostatic deflector. 



II. EXPERIMENTAL SETUP 

A. Molecular beam apparatus 

An exploded schematic of the molecular beam appara- 

It consists 



tus, used in these studies, is shown in Fig. 1 



of three vacuum chambers: the source chamber, pumped 
by a 2000 1/s turbomolecular pump, followed by the de- 
flector and the target chamber which are both pumped 
by 500 1/s turbomolecular pumps. The molecular beam 
is formed by expanding a mixture of carrier gas (typi- 
cally 90 bar of Helium) and the molecular gas (~5 mbar 
of benzonitrile) into the source chamber through a high- 
pressure Even Lavie valve [26l |27] heated to 40 °C. Op- 
erating the valve at high stagnation pressure and fine 
tuning the valve opening time ensures formation of a 



B. Electrostatic deflector 

The deflector consists of a trough-shaped electrode 
with a 3.2-mm-inner-radius of curvature kept at ground 
potential and a polished rod with a radius of 3.0 mm 
to which high voltage (up to 10 kV) can be applied [29]. 
The rod is kept floating above the trough by two pieces of 
macor, providing a minimum distance of 0.9 mm between 
the electrodes and 1.4-mm-separation at the molecular 
beam axis. This electrode geometry creates a two wire 
field, that is strongly inhomogeneous along the (vertical) 
Y-axis while it is almost homogeneous in the X direc- 
tion [30]. Therefore, polar molecules traversing the de- 
flector experience a force in the vertical direction with 
molecules in high-field seeking states being deflected up- 



wards. The extent of deflection for a molecule depends 
on the effective dipole moment of the particular molec- 
ular quantum state. As the molecules exit the deflec- 
tor, they enter a parallel plate capacitor consisting of 
two 17-cm-long polished stainless steel plates separated 
by 2.7 mm. The plates are oriented parallel to the de- 
flection direction. The capacitor creates a constant dc 
electric field on the order of 2 kV/cm along the path 
of the molecules moving from the deflector to the laser 
interaction zone. This finite field prohibits population 
transfer between quantum states due to Majorana tran- 
sitions between M-levels of individual rotational states 
at zero field, and it mitigates transfer due to diabatic 
traversal of avoided crossings induced by rotations or in- 
homogeneities of the electric field, which occur frequently 
at small field strengths [3TH33] 



50 mm-diameter 30 cm- focal- length spherical lens. The 
probe pulse is electronically synchronized to the peak of 
the YAG pulse. 

Using telescopes in each laser beam the spotsizes are 
adjusted so that the probe beam is more tightly focused 
than the YAG beam {ouj^^ = 34 /im and cjP'°^^ = 
21 /im, measured by scanning a 10 /im pinhole). The 
intensity of the probe pulse is adjusted by inserting thin 
neutral density filters in the beam, the intensity of the 
YAG pulse by rotating a half wave plate placed between 
two thin-film polarizers, upstream, and one thin film 
polarizer downstream. The polarization state of either 
beam is controlled by a half- wave and a quarter- wave 
plate inserted in each beam just before the final degree 
overlap mirror. 



C. Velocity map imaging spectrometer and laser 
beams 



III. EXPERIMENTAL RESULTS 



A. Deflection of benzonitrile 



The target chamber houses a velocity map imaging 
(VMI) spectrometer consisting of an open three-electrode 
electrostatic lens \M . Between the first two electrodes 
(repeller and extractor) the molecular beam is crossed 
by two (or more) pulsed laser beams as indicated with 
the red and blue (and green) cylinders in Fig. 1 Apply- 
ing positive (negative) voltages to the electrostatic lens 
ions (electrons) in the interaction region are accelerated 
towards a 50-mm-diameter microchannel plate (MCP) 
backed by a similar size phosphor screen. The ion (or 
electron) images are recorded by a CCD camera monitor- 
ing the fluorescent phosphor screen, and on-line software 
analysis determines and saves the coordinates of each in- 
dividual particle hit. When imaging ions the MCP de- 
tector can be gated by applying a fast high- voltage pulse 
to its front thereby obtaining images for a selected range 
of mass-to-charge ratios. The spectrometer is caged in a 
single concentric mu-metal cylinder to minimize the in- 
fluence of external magnetic fields on the trajectories of 
the electrons. 

The laser systems, as well as the spatial and temporal 
overlap of different laser beams with the molecular beam 
is similar to that reported before [28 so the description 
here is brief. The fundamental output (1064 nm) from 
a 20 Hz, NdiYAG laser is used to adiabatically align the 
molecules. The pulse duration is 10 ns. Part of the out- 
put from a 1 kHz amplified Tiisapphire laser (800 nm) is 
compressed to 30 fs and used to either Coulomb explode 
the molecules [See. HfB^ , for detecting their spatial orien- 
tation or to singly ionize the molecules for creating photo- 
electrons ( Sec. HIC]) or f or characterizing the molecular 
beam profile (Sec. HI A). The alignment beam (here- 
after termed the YAG beam or the YAG pulse) and the 
fs beam (hereafter termed the probe beam or the probe 
pulse) are collinearly overlapped using a zero-degree mir- 
ror (high reflectivity at 1064 nm and high transmission 
at 800 nm) and focused onto the molecular beam by a 



The effect of the deflector on the molecular beam is il- 



lustrated by the vertical intensity profiles in |Fig. 2| These 
profiles are obtained by recording the BN+ signal [BN = 
benzonitrile (C7H5N)] from photoionization due to the 
femtosecond probe laser (/probe ^ 10^^ W/cm^) as a 
function of the vertical position of the laser focus. When 
the deflector is turned off, the molecular beam extends 
over ~ 1.5 mm, mainly determined by the diameter of 
the skimmer before the deflector and the end-aperture of 
the deflector. When the deflector is turned on, the molec- 
ular beam profile broadens and shifts upwards (positive 
Y values). Molecules in the lowest rotational quantum 
states have the largest Stark shifts and are, therefore, de- 
flected the most [21, 29 . In the measurements described 
below, experiments were conducted on these quantum- 
state selected molecules simply by positioning the laser 
foci close to the upper cut-off region in the 5 k V profil e 
(Y = 1.55 mm), as indicated by the red arrow in Fig. 2 



B. Alignment and orientation 

The molecules studied here are asymmetric tops. Their 
alignment and orientation is controlled by the combina- 
tion of the YAG pulse and the weak static electric field in 
the VMI spectrometer. Based on previous experiments 
and theory a linearly polarized YAG pulse is used to in- 
duce ID alignment and orientation whereas 3D alignment 
and orientation is created by applying an elliptically po- 
larized YAG pulse [H EH ESHSg. Here ID alignment 
refers to confinement of a single molecular axis along the 
YAG polarization axis. Due to the nature of the polar- 
izability interaction between the molecule and the lin- 
early polarized alignment field it is the most polarizable 
axis that is aligned [14 . For benzonitrile, discussed here, 
this is the C2-axis (i.e., the C-N axis). Similarly, 3D 
alignment refers to confinement of three perpendicular 




(a) Repeller Extractor Ground 



0.0 0.5 

Y(mm) 

FIG. 2: (Color online) Vertical profiles of the molecular beam 
measured by recording the BN^ signal as a function of the 
vertical position of the fs probe beam focus. The experimental 
data are shown by black squares (deflector off, kV) and red 
circles (5 kV). The connecting lines serve to guide the eye. 
The arrows indicate the position of the laser foci for acquiring 
ion and electron images of undeflected (black), Y = mm, and 
deflected (red), Y = 1.55 mm, molecules. 



molecular axes [IT. In practice, the elliptically polarized 
alignment field confines the most polarizable axis along 
the major polarization axis, the second most polarizable 
axis (for BN the axis in the molecular plane perpendic- 
ular to the C2-axis) along the minor polarization axis 
and, consequently, the least polarizable axis (for BN the 
axis perpendicular to the molecular plane) perpendicu- 
lar to the polarization plane. Finally, orientation refers 
to whether the dipole moment of the molecules point in 
a particular direction with respect to a laboratory-fixed 
reference direction, which here is the weak static electric 
field in the VMI spectrometer. 

To test that alignment and orientation of BN occurs, 
the molecules are Coulomb exploded by the intense fs 
probe pulse (/probe = 5.4 x 10^^ W/cm^) and recoil- 
ing ionic fragments are recorded by the 2D ion detector, 
a method used previously for many other molecules in 
our laboratory. For a planar molecule of C2v symme- 
try, like BN, our standard procedure is to first estab- 
lish that the C2-axis is aligned. This is accomplished by 
recording C N+ ion s because they primarily recoil along 
the C2-axis. Fig. 3 ^b) shows a CN+ image recorded with 
the probe pulse only, linearly polarized in the detector 
plane. A mild confinement of the ions is observed re- 
flecting the enhanced ionization probability for molecules 
that initially have their C2-axis pointing along the probe 
pulse polarization. The confinement is quantified by 



(cos ^2d) =0.63, where O2B is the angle between the pro- 
jection of the CN+ recoil velocity on the detector plane 
and the YAG polarization [28]. When the YAG pulse is 
included, linearly polarized vertically, corresponding to 
P = 90° {P denoting the angle between the YAG polariza- 
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FIG. 3: (Color online) (a) Schematic of the velocity map 
imaging spectrometer used to detect ions or electrons. The ID 
alignment of the BN molecules is determined by the alignment 
laser polarization, here shown for /3 — 45°. The static electric 
field of the spectrometer, pointing from the repeller to the 
extractor electrode for ion detection, breaks the head-for-tail 
symmetry by preferentially placing the CN-end towards the 
repeller where the positive voltage is highest. When detect- 
ing electrons the polarity of the electrodes is inverted forcing 
the CN-end of the molecules towards the extractor electrode. 
(b)-(f) CN^ ion images, (b): probe pulse only, vertically po- 
larized, (c), (d) P — 90°. In (d) the rotational temperature 
is significantly higher than in (c) and no quantum state selec- 
tion occurs (see text), (e), (f): 13 — 45° and -45°, respectively. 
For all images Fstat = 286 V/cm . In (e) / (f) 71 % / 28 % 
of all CN^ ions appear in the upper half of the detector. 



tion and the static e lectric field), the angular confinement 

and (cos^ 6>2d) increases to 



sharpens markedly Fig. 3 



0.89. This shows that the BN molecules are ID aligned 
with their C2-axis along the YAG polarization. The 
(cos^ ^2d) value is quite a bit smaller than that observed 
for ID alignment of iodobenzene ((cos^^2d) = 0.97) 
[2T] under similar conditions of initial rotational cool- 
ing, quantum-state selection and alignment pulse inten- 
sity. This may seem surprising since the components of 
the polarizability tensor for BN are similar to those of 
iodobenzene and thus the two molecules should exhibit 
approximately the same degree of ID adiabatic align- 
ment. The reason for the lower (cos^^2d) value for BN 
is most likely that CN+ ions are not ideal observables 
for deducing the spatial orientation of the C2-axis be- 
cause, unlike 1+ ions used in the iodobenzene characteri- 
zation, it overlaps with other fragment ions in the time of 



flight spectrum of C7H5N. Hence, the ion images of |Fig. 3 



also contain contributions from hydrocarbon fragments 
from the benzene ring of equal (or nearly equal) mass- 
to-charge-ratio (e.g. C2H2"). These "contaminant" ions 
are not expected to recoil along the C-CN axis. In the 
alignment measurements we try to minimize the contri- 
bution from the contaminant ions by only analyzing the 
CN+ ions with the highest kinetic energy corresponding 
to the outermost ring [39] seen in the images of Fig. 3 



However, just a small contributions from the contami- 
nant ions ejected in directions away from the YAG po- 
larization is enough to signiflcantly lower the observed 
(cos^ ^2d) value. 



For reasons discussed in Sec. |HIC| alignment was also 
measured for conditions of higher rotational temperature 
and no rotational state selection. In practice, this is 
achieved by lowering the He stagnation pressure from 85 
to 15 bar [28 and, in addition, turning off the deflector 
voltages. The pronounced decrease in the degree of align- 
ment is clear by comparing Fig. pfd) (15 bar He, deflector 
turned off) with Fig. [3](c) (90 bar He, selection of lowest 
rotational quantum states). Quantitatively, (cos^^2d) is 
0.76 and 0.89 in Figs. [31(d) and (c), respectively. 

To achieve 3D alignment the polarization of the YAG 
pulse is changed from linear to elliptical. When an el- 
lipticity of 3:1 is chosen, the value referring to the ratio 
of the intensities measured along the major and the mi- 
nor polarization axis, the CN+ ion image (not shown) 
recorded with the major polarization axis in the plane 
of the detector (vertical) is very similar to the image in 



Fig. 3| (c), showing an equally strong degree of ID align- 
ment of the C2-axis. To additionally conflrm conflne- 
ment of the molecular plane H+ ions were recorded, see 
|Fig. 4[ following an approach employed in previous stud- 
ies [38^, "W of 3D alignment. When the probe laser only 
is used the image is circularly symmetric [Fig. [4] (a)]. In- 
cluding the linearly polarized YAG pulse, polarized per- 
pendicular to the detector plane, does not change the 
circular symmetry [Fig. p](b)]. This is consistent with 
the C2-axis being aligned along the YAG polarization 
vector and the molecular plane free to rotate around the 
arrested axis. When the 3:1 elliptically polarized YAG 
pulse, its major polarizability axis perpendicular to the 
detector, is used, the circular symmetry is broken and, 
instead, the H+ ions localize around the vertical axis cor- 
responding to the minor axis of the polarization ellipse. 
Similar to previous studies [38l HO] we interpret this as 
conflnement of the molecular plane with the second most 
polarizable axis (perpendicular to the C2-axis but still in 
the molecular plane) aligned along the minor axis of the 
ellipse. 

To show orientation of BN the polarization of the lin- 
early polarized YAG pulse is rotated to /3 = 45° (—45°) 
meaning that the C2-axis of the molecules is positioned 
at 45 (-45) with respect to the static_electric fleld, Fgtat, 

a)]. The resulting 
Sle) and (f). The 



of the VMI spectrometer [see Fig. |3 

CN+ images are displayed in Figs. 

appearance of more CN+ ions on the lower (upper) part 




FIG. 4: Ion images of H^ recorded with the probe laser 
(5.4 X 10^^ W/cm^) linearly polarized perpendicular to the 
detector plane, (a) No YAG pulse, (b) YAG pulse included, 
linearly polarized perpendicular to the detector, (c) YAG 
pulse included, elliptically polarized. The major axis is per- 
pendicular to the detector and the ellipticity-intensity ratio 
is 3:1. The center has been removed to suppress a large spike 
of very low kinetic energy H^ ions stemming from ionization 
of residual molecules in the target chamber. 



of the detector for 45 (-45) shows that the molecules are 
oriented with the CN-end preferentially pointing towards 
the repeller electrode. Hence, the permanent dipole mo- 
ment of BN is pointing in the direction from the repeller 
to the extractor as expected. The up-down asymmetry 
remains unchanged when an elliptically rather than a 
linearly polarized YAG pulse is employed to align the 
molecules. These flndings of alignment and orientation 
are fully consistent with recent studies on iodobenzene 
[211 [29] and 2,6-difluoroiodobenzene p^. 

The degree of orientation is characterized by the ratio 
of ions detected on the upper part to the total number of 
ions detected. For /3 = 45° it is 27 % and for /3 = -45° it 
is 71 %. For several reasons these numbers underestimate 
the degree of orientation utilized in the photoelectron 
experiments described below. The presence of contam- 
inant ions with the same mass-to-charge ratio as CN+ 
will reduce the up-down asymmetry. In addition, the ion 
images are recorded at Fstat=286 V/cm which is signifl- 
cantly lower than the electrostatic extraction fleld in the 
photoelectro n angular d istribution (PAD) measurements 
presented in Sec. IIIC| (Fgtat = 467 V/cm). In fact, the 
increase of the static fleld will be even larger because the 
ion measurements are recorded at /3 = 45° , i.e., the ef- 
fective value of Fstat is reduced by >/2. By contrast for 
the PAD measurements the alignment fleld is polarized 
parallel to Fgtat- As a consequence, we estimate that 
the orientation in the PAD geometry corresponds to at 
least 80 % of the molecules have their CN-end facing the 
extractor electrode. 



C. PADs from single ionization of benzonitrile 



For the PAD experiments the same experimental setup 



as for the ion detection described in Sec. IIIB is used. 



but some essential parameters are changed. The po- 
larization state of the 30 fs probe pulses is changed 
from linear to circular and the intensity is lowered to 
::^ 1.2 X 10^^ W/cm^ corresponding to a regime where 



the BN molecules only undergo single ionization with 
essentially no fragmentation. Also, the polarization of 
the alignment pulse is changed such that the major axis 
is parallel to the static field axis. Furthermore, to ex- 
tract electrons instead of ions in the PAD measurements 
the polarity of the velocity map imaging spectrometer is 
inverted. Hereby, the BN molecules are confined along 
the static field axis with the CN-end facing the extractor 
electrode. 

The photoelectron images are shown in |Fig. 5[ Ap- 
plying only the probe pulse [Figs. [sJa.S) and (a. 4)] the 
electrons emerge in a stripe parallel to the polarization 
plane (Y,Z) of the probe pulse for both left and right 
circularly polarized (LCP and RCP) pulses. A weak up- 
down asymmetry is observed, that reverses as the he- 
licity of the probe pulse is fiipped. The asymmetry is 
more easily seen in the angular distributions [Figs. [5|a.2) 
and (a. 5)], obtained by radially integrating the respective 
images from momentum 0.02 to 0.8 a.u. When the YAG 
pulse, linearly polarized along Fgtat, is included to induce 
strong ID alignment and orientation, schematically illus- 
trated in [Fig. 5| b.l), the up-down asymmetry increases 
[Figs. [sJb.S) and (b.4)]. For LCP (RCP) probe pulses 
the number of electrons detected in the upper part com- 
pared to the total number in the image is ~55% (44%). 
The enhanced asymmetry is also clearly visible in the 
angular distributions [Figs. |5|b.2) and (b.3)]. As ex- 
plained in |Sec. IV| the up-down asymmetry results from 
the orientational dependence of the ionization yield and 
the force on the electrons, by the probe laser field, fol- 
lowing their detachment from the molecule. In particu- 
lar, we show that ionization of molecules, ID aligned and 
oriented with their dipole moment along the direction of 
Fstat 1 leads to electrons with an upward (downward) mo- 
mentum for LCP (RCP) probe pulses in agreement with 
the observations. The observation of an asymmetry with- 
out the YAG pulse results from mild orientation due to 
the interaction between the permanent dipole moment 
and the weak static field and/or state selection by the 
deflector [38 . 

When the BN molecules are 3D aligned and oriented, 
obtained by using an elliptic ally ra ther than a linearly po- 
larized YAG pulse (see Se cJIIB ), striking substructures 
appear in the PADs [Figs. |5|c.3) and (c.4)]. Notably, the 
electron emission is suppressed along the vertical direc- 
tion, clearly seen as a dip in the corresponding angular 
distributions [Figs. |5|c.2) and (c.5)] at a = 90° and 
270°, where a is the angle measured counterclockwise 
with respect to the X direction [see Fig. [5|a.3)]. The an- 
gular distributions also show that the maximum in the 
electron emission occurs at an angle 



for the 3D aligned and oriented molecules, corresponds 
to the (Y,Z) plane. Regarding the up-down asymmetry 
of the electron distributions [panels (c.2) and (c.5)] it 
is essentially the same as for the electron distributions 
obtained for the ID aligned and oriented molecules [Figs. 
(b.2) and (b.5)]. This is in accordance with expectations 
since the degree of ID alignment and orientation should 
be essentially the same for the linearly and the elliptically 
polarized YAG pulse (verified experimentally). 

To test that the observed structures in the PADs are 
not due to perturbation from the YAG pulse on the ion- 
ization process by the fs probe pulse, experiments were 
carried out on molecules with a much lower degree of 
alignment but identical YAG pulse intensity. This is 
achieved by employing molecules with higher rotational 



l]p 



18° 



(1) 



away from vertical. As explained in Sec. IV the local 
minimum of electrons along the vertical direction results 
from suppression of electron emission in the nodal plane 
of the highest-occupied molecular orbital (HOMO) [and 
the first orbital below HOMO in energy, HOMO-1] which. 



temperature, achieved as explained in Sec. IIIB Figure 
|6] shows the photoelectron images obtained for LCP and 
RCP when the molecules are 3D aligned and oriented. 
Images similar to those in Fig. [5] are kept as a reference 
[Figs. [6Fa.2) and (b.2)]. Both the up-down asymme- 
try and the lobe-like structures are completely absent in 
Figs. |6ja.l) and (b.l), thereby showing the importance 
of a tightly aligned and oriented target for observation of 
the effects in the PADs. The difference between the ex- 
periments on strongly versus weakly aligned and oriented 
molecules is also obvious from the angular distribution 
displayed in Figs. |6Fa.3) and (b.3). 



IV. THEORY 

A. Stark-shift corrected tunneling theory 

To describe the response of molecules with large po- 
larizabilities and dipole moments to strong, near infrared 
laser pulses requires the Stark shifts of the energy levels 
be included, see Ref. [23land the Supplementary Informa- 
tion in Ref. |22l In Ref. |23l we modified the tunneling 
model to include the Stark shifts of the molecule and its 
cation at the instantaneous value of the external field. 
We briefiy summarize the model here for completeness. 

When the Stark shifts are taken into account, the field- 
free ionization potential /p(0) is modified and reads (up 
to second order in the field strength F) 



/^(F) = /^(O) + Am • F + ^F^AaF, 



where 



A/L6 = ^ 



M 



M 



Aa = OL 



M 



OL 



(2) 



(3) 



II is the permanent dipole moment, a is the polarizabil- 
ity tensor, and the superscripts M and / are referring to 
the molecule and the unrelaxed cation, respectively. We 
use the unrelaxed cation, since the nuclei have no time to 
move to the new equlibrium position during the tunneling 
process. At the field strengths of concern the shifts in- 
duced by the hyperpolarizability and higher-order terms 
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FIG. 5: Two-dimensional momentum images of electrons produced when BN molecules are ionized by a LCP probe pulse or 
a RCP probe pulse. The corresponding angular distributions in the detector plane are shown next to the images. In row (a) 
the molecules are essentially randomly oriented (no alignment pulse). In row (b) the molecules are ID aligned and oriented 
(linearly polarized alignment pulse) and in row (c) the molecules are 3D aligned and oriented (elliptically polarized alignment 
pulse). The intensity of the probe pulse is 1.2 x lO^'^ W/cm^. The intensity of the YAG pulse is 7 x 10^^ W/cm^. In row (c) 
the intensity ratio between the major and minor axis is 3:1. The (X,Y,Z) coordinate system is oriented identical to the one 
shown in Fig. ^ 
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FIG. 6: Two-dimensional momentum images of electrons produced when BN molecules are ionized by a RCP probe pulse (a.l), 
(a.2) or a LCP probe pulse (b.l), (b.2). The images in column 1 (2) are recorded for weakly (strongly) 3D aligned and oriented 
molecules obtained by using 15 (85) bar of He stagnation pressure (see text). The laser pulse parameters were the same as 
those used in Fig. Is] The angular distributions in the detector plane are compared in (a. 3) and (b.3). 



are negligible compared to the linear and quadratic terms 
accounted for in Eq. ([2|. We note that the dipole mo- 
ment and polarizability differences in Eq. ^ are exactly 
the dipole and polarizability of the orbital actively in- 
volved in the tunneling process, e.g., the HOMO and 
HOMO-1 in the present case. Hence the discussion could 
equivalently be performed referring to the orbital dipole 
moment and polarizability. Such an approach is natu- 
ral in the case of calculations with the time-dependent 
Schrodinger equation within the single-active electron ap- 
proximation [41,. 

In the tunneling regime, i.e., when the Keldysh pa- 
rameter 7 [42] is smaller than unity, the static Stark 
shift should be accounted for and the standard tunnel- 
ing model of Ref. [43 modified accordingly by using the 
shifted ionization potential of Eq. Q. These changes 
are crucial in order to reproduce the experimental asym- 
metry in the photoelectron momentum distribution [22 . 
First, for simplicity we restrict for simplicity our con- 
sideration to the case were the target molecule can be 
modeled by an s-orbital. We include a factor to account 
for the saturation of ionization at high intensities [44], 
and the Stark-corrected tunnel ionization rate (up to a 
multiplicative constant) becomes 



w{F) 




2K{Fy 
3F 



where 



/.(F) 



2/JF) 



(4) 



(5) 



We note from Eqs. ([2| and ([5| that k,(F) is a function of 
the strength and the orientation of the field with respect 
to the molecular frame. 

The Keldysh parameter for the current experiment in- 
volving benzonitrile amounts to 1.17. It has been shown 
that for increasing molecular size the ionization becomes 
tunnel-like at larger Keldysh parameters [45 ; therefore 
we use the tunneling model here. In benzonitrile the 
tunnel exit is far from the centre of mass, so in circularly 
polarized laser fields, where no rescattering occurs, the 
influence of the molecular potential on the escaping elec- 
tron can be neglected. Accordingly, after the initial ion- 
ization step, the subsequent propagation of the electron 
can be estimated accurately by solving the classical equa- 
tions of motion in the laser field only [46 . We consider 
a left circularly polarized (LCP) probe pulse, containing 
more than 10 cycles, which is switched off adiabatically 
at large times. These assumptions allow us to consider a 
field with constant amplitude 



F{t) = Fq {sin{ujt)ey + cos(cjt)ez) . 



(6) 



Here Fq is the peak field strength and uj the angular 
frequency. The final photoelectron momentum at the 2D 
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FIG. 7: The HOMO and HOMO-1 of benzonitrile (C7H5N). 
The molecular plane coincides with the (Y,Z) plane, and the 
detector is on the right of these images. The insets in dashed 
lines, show the view of the orbit als from the detector, show- 
ing that both orbitals have nodal planes in the polarization 
plane of the circularly polarized laser pulses. The orbitals and 
ionization potentials have been obtained by quantum chem- 
istry caculations in GAMESS 47 at the Hartree-Fock level 
of theory in conjunction with the valence triple-zeta basis set. 



detector screen can now be straightforwardly related to 
the vector potential at the instant of ionization to, that 
is 



Py =- Ayiujto) 
Pz =- Az{ujto) 



u 



■ cos(cjto) 



— sin(cjto) 



(7) 



where ujto is the angle between the Z-axis and the field 
vector at the instant of ionization. Using Eqs. Q and 
(JTI, for each field orientation, i.e., to, the final momenta 
in the (Y,Z) plane and the associated tunneling proba- 
bility can be calculated. 



B. Modeling the off— the-nodal plane photoelectron 
emission by tunneling theory 

The structure of the molecular orbitals of benzonitrile 
with the lowest binding energies, combined with the par- 
ticular 3D orientation in the present experimental setup, 
result in the unique X-shaped pattern in the electron 
momentum distribution shown in Fig. [5J and the off- 
the-nodal plane emission angle reported in Eq. ([l]). To 
explain such a pattern requires a more advanced treat- 
ment than the one presented in Ref. 22 , The HOMO and 
HOMO-1 of benzonitrile are shown in Fig. [7[ The most 
important feature for the present development is that the 
relevant experimental orientation of the molecule is such 
that the nodal planes of both the HOMO and HOMO-1 
coincide with the molecular plane which in turn coin- 
cides with the polarization plane of the laser pulse. In 
other words, neither HOMO nor HOMO-1 have electron 
density in the polarization plane of the pulse ([6|. This 



geometry poses a theoretical challenge in tunneling mod- 
els, since conventional tunneling and molecular tunneling 
theory is inapplicable. The simple reason for their fail- 
ure is that the standard assumption, that tunneling of 
the negatively charged electron happens in a very narrow 
cone around the direction opposite to the field [43l l48l - 
[50], is no longer valid. In that direction there is no bound 
state wave function to overlap with the continuum emerg- 
ing from the outer turning point, and no ionization will 
occur due to the nodal plane. In short, in the nodal 
plane there is no population which can tunnel, and conse- 
quently emission occur off-the-nodal plane. The simplest 
cases where this is evident is tunnel emission from a hy- 
drogenic 2p^ or 2px orbital when the field is aligned along 
the Z-axis. These states are linear combinations of the 
states |n = 2, / = 1, 772 = 1) and |n = 2, / = l,77i = —1). 
According to tunneling theory, the contributions from 
states with m = 1 and m = — 1 should be added incoher- 
ently. Then, following Ref. [51], the transverse momen- 
tum distribution is proportional to 



exp 




(8) 



and therefore from such a state one would observe a peak 
at px = in the momentum distributions recorded on 
the detector screen in this experiment. However, as is ev- 
ident from the present experiment, the emission occurs 
off-the-nodal plane. Inclusion of off-the-nodal plane emis- 
sion in the tunneling model is nontrivial, since the dom- 
inant photoelectron momenta under tunneling emission 
(momentum distribution of tunneled electron at its birth) 
should now depend on the particular initial state. For ex- 
ample, 2p^ or 2px orbit als would have different transverse 
momentum distributions. Concerning the derivation of 
the tunneling rate, we can, in this situation, no longer 
separate the azimuthal angle around F(t) and thus can- 
not add the individual m components incoherently. We 
note here that even when static Stark shifts are included, 
the information about the nodal plane of the initial state 
survives and cannot be washed out by the polarization 
effects. 

Motivated by earlier work on tunneling theory for large 
systems [52 , and by the fact that the saddle-point for- 
mulation of the length-gauge strong-field approximation 
(SFA) [42 in the limit of small Keldysh parameter agrees 
very well with the tunneling theory [53 (at least for 
a zero-range potential), we modify the transverse mo- 
mentum distribution of Eq. ([8| by multiplying it by the 
squared Fourier transform of the initial wavefunction. 
More precisely, for the 2px orbital that means multipli- 
cation of the transverse momentum distribution ([8| by 
the factor (ppcos(0))^, i.e.. 



^(jPpA) --PpCos^((/>)exp 






(9) 



where Pp and (j) are the 2-dimensional polar coordinates 
in the (X,Y) plane. In Eq. ([9| (ppcos(0))^ is the dom- 



inating factor relating to the angular behavior of the 
squared Fourier transform of 2px at p^ = 0, and a factor 
{{p/ n^ + 1)~^ has been neglected. Finally, in the above 
equation for the transverse momentum distribution, we 
have used the field-free ionization potential (appearing 
as a square root in the exponential) since, as we have 
checked, adding Stark shifts to this part of the distribu- 
tion does not change the results much. 

Using the distribution ([9| the expectation value of the 
X-component of the transversal momentum (we assume 
the field is oriented along the Z-axis), px = PpCos((/)) in 
the interval px > can be calculated analytically as 



p7r/2 



{p. 



'X 



J^ dppPl 111^2 ^^ cos(0)^(pp, (j)) 
ir dPpPp Jli% d(t)w(jpp, 0) 




(10) 

where n = ^J2Lp(fi). We note that the expectation value 
of py for the tunneled electron at birth is {py) = and 
that the probability of emission with px = is zero, i. e., 
there is no emission in the nodal plane of the initial state 
although the field is along the Z-axis. The final momenta 
of the tunneled electron along the Y- and Z-axis are de- 
fined by Eq. ([tI and it follows that the maximal final elec- 
tron momentum along the Y-direction is p'^^^ = Fq/uj. 
We can define an angle Vt = arctan((j9x)/Py'^^) as a mea- 
sure for the off-the-nodal plane emission. With the above 
expressions, this angle can be calculated analytically and 
the results reads 

Q = arctan((px)/Py'''^) = arctan ( . 



arctan | —^^y 




(11) 



where 7 = uok/Fq is the Keldysh parameter [42]. The 
above expression is, to our knowledge, the first quan- 
tification of the angle of off-the-nodal plane emission as 
a function of the field intensity and frequency and the 
binding energy of the atomic/molecular system. It has 
the intuitively correct behavior with respect to increase of 
the field intensity and the binding energy: the angle de- 
creases. Since both peak field Fq and n appear as square 
roots in the expression, in the limit of large intensities, 
where also the tunnel model is more valid, the angle is 
inversely proportional to the fourth root of intensity, i. e., 
its dependence on intensity is very weak. 

To test our model, we imagine a fictitious experiment 
where we have an 'oriented' atom in a p^^ orbital and the 
field- free ionization potential /p(0) = 0.5 a.u.. We as- 
sume that, identically to the experiment involving ben- 
zonitrile, ionization is induced by a LCP field in the (Y,Z) 
plane and that the momentum distribution is recorded on 
a detector coinciding with the (X,Y) plane. The final mo- 
mentum distribution recorded on the screen can be cal- 
culated for this case using the above off-the-nodal-plane 
tunnel emission model. We have compared the results 
from our model to the calculations within the strong- 
field approximation [42] assuming a 10 cycle circularly 
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FIG. 8: Expectation value of (a) px and (b) off-nodal-plane 
emission angle Q as a function of intensity of the circularly 
polarized field for a px orbital with binding energy of 0.5 a.u. 
similar to the HOMO of benzonitrile and a laser wavelength 
of 800 nm. The full curve denotes the present analytic the- 
ory, and the dots denote the results obtained from calculation 
within the strong- field approximation [42, . 



polarized laser pulse with a sine-square envelope, defined 
in terms of the vector potential A(t) interacting with 
the same 2px state. The comparison between the pre- 
dictions of (px) and Q by the two methods is given in 
Fig. [8J It can be seen that for large intensities, the tun- 
neling model and the SFA calculations largely agree for 
both quantities. As the intensity of the field decreases, 
or equivalent ly, the Keldysh parameter 7 increases, the 
discrepancy between the SFA and the tunneling model 
results increases, reflecting the gradual departure from 
the tunneling regime. 



V. STRONG-FIELD IONIZATION OF 
3D-ORIENTED BENZONITRILE 

Armed with our generalized tunneling model including 
both the static Stark shifts and the off-the-nodal-plane 
emission, we now address the experimental results for 
benzonitrile. 

As illustrated in Fig. [5| two sets of experiments 
have been performed involving pre-oriented benzonitrile 
molecules: one with ID-oriented molecules and the other 
with 3D-oriented molecules. The results obtained with 
ID and 3D orientation show almost the same degree of 
up-down asymmetry (Fig. [5|. The asymmetry can be 
ascribed to the dependence of the ionization pote ntial o n 
the static Stark shifts of the energy levels (see Sec IV A), 
i.e., the small difference between the dipole moment of 
benzonitrile and its cation (Appendix A). The case of ID- 
oriented benzonitrile can be treated theoretically by the 
same method as introduced in Ref. 22 for the description 
of the linear DCS molecule. In the ID oriented case, the 
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FIG. 9: (a) Momentum distributions including volume effects 
for benzonitrile at 1.23 x 10^^ W/cm^ for ionization from the 
HOMO and HOMO-1. (b) Angular distributions derived from 
momentum distributions. Black (full): HOMO+HOMO-1, 
red (dashed): HOMO, blue (dotted): HOMO-1. The results 
should be compared to the corresponding experimental find- 
ings in Fig. [5] 



molecules are allowed to rotate freely around its molecu- 
lar axis and the effect of off-the-nodal-plane emission is 
washed out. In the experiment on 3D aligned and ori- 
ented benzonitrile, on the other hand, the polarization 
plane of the circularly polarized laser lies entirely in the 
nodal plane of both HOMO and HOMO-1. Indeed, in 
this case, off-the-nodal-plane emission occurs and mani- 
fests itself as an X-like pattern in the photoelectron mo- 
mentum distribution on the detector screen, and here 
we focus on the 3D-oriented benzonitrile. To describe 
the experiment, we use the theory for off-the- nodal plane 
emission, developed in the previous section. The polariz- 
ability of benzonitrile is very large (see Apendix A) , so we 
expect that the initial orbital would be modified by the 
polarization response of the system, and we assume that 
any structure of the orbital off the nodal plane will be 
washed out. Thus, when the Stark shift is not included, 
the total probability of tunnel emission will not depend 
on the (finite) angle of the electric field vector with re- 
spect to the Z-axis. Therefore, we model all orbitals of 
benzonitrile that have their nodal plane coinciding with 
the polarization plane with the simplest possible orbital 
having such a nodal plane: the p^^ orbital. All theoretical 
results presented here assume perfect 3D orientation of 
the molecule. Using Eqs. Q, ([7| and ([9|, the momentum 
distributions at a particular intensity can be calculated. 
We perform focal volume integration by including inten- 
sities ranging from 1.2 x 10^^ W/cm^ to 0.5 x 10^^ W/cm^ 
with a step of 5 x 10^^ W/cm^. The resulting momentum 
distributions from HOMO, HOMO-1 and the combined 
response of HOMO and HOMO-1 are shown in Fig. |9] 

As seen by comparison of Fig. J5] and F ig. |9) at larger 
radial final momenta Pp = \/v\^Vy-> there is good 
agreement between the experimental and theoretical dif- 
ferential momentum distributions. Also, the agreement 
between theory and experiment for the angle off-the- 
nodal-plane emission ^ is good: the experiment yields 
- 18 degrees (Fig, [sj [Eq. ([l])]) and the theory 18.8 de- 
grees [Eq. (11)]. In agreement with the calculation the 
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suppression of electron emission occurs at a = 90 and 270 
degrees. The degree of suppression seen experimentally 
is, however, much less pronounced that the complete sup- 
pression, i.e. no electron emission with px = 0, predicted 
theoretically. One factor contributing to this disagree- 
ment is the lack of perfect alignment of the molecular 
plane in the experiment. Another source that can pro- 
duce a signal at px = is emission from the HOMO-2 
that is bound with additional 0.1 a.u. compared to the 
HOMO. The HOMO-2 has population in the polariza- 
tion plane of the field so that its contribution can still be 
significant (see Fig. 4 in the Supplementary Information 
ofRef. [22J). 

The present theory is capable of describing the small 
up/total asymmetry observed in the experiment by at- 
tributing it to the Stark shift of the ionization potential 
due to the dipole and polarizability of the HOMO orbital, 
i.e., the difference of the permanent dipole moment of the 
molecule with respect to the permanent dipole moment 
of the unrelaxed cation. The calculated up/total asym- 
metry for the case of 3D-oriented molecule at 1.23 x 10^^ 
W/cm^ is 0.52 in fair agreement with the experimental 
results. 



VI. CONCLUSIONS 

We have studied molecular frame photoelectron distri- 
butions (MFPADs) following strong field single ioniza- 
tion, with circularly polarized, nonresonant fs pulses, of 
benzonitrile molecules aligned and oriented in space by 
mixed laser and static electric fields. The most striking 
experimental observation is suppression of electron emis- 
sion in directions that correspond to nodal planes of the 
orbitals from which ionization occurs. The experimental 
findings are explained using a modified tunneling the- 
ory. Our theoretical analysis shows that the circularly 
polarized ionizing field drives the photoelectrons away 
from the cation, minimizing effects of post-ionization in- 
teraction, and leads to a clear signature of the orbital 
structure of the molecule at the time of ionization. In 
particular, we have shown here that nodal structures can 
be mapped out, a possibility recently examined further 
theoretically [54 . 

An important future challenge is to extend the present 
work to the study of time-dependent phenomena. This 
may be realized, for instance, by a pump-probe scheme 
where the probe pulse promotes molecules from the elec- 
tronic ground state to an electronically excited state. 
One class of such photoinduced reactions that could be 
well suited for the strong field probing, introduced here, 
concerns large amplitude vibrational motion in the ex- 



cited state. This occurs for, instance, in axially chiral 
systems such as biphenyl or phenyl-pyrrole [55 where 
the dynamics in the excited states corresponds to tor- 
sional motion of two aromatic rings. Alignment as well 
as orientation of such systems is feasible j56l IS?] and 
the present probe method could map out the changing 
valence electronic structure during torsional motion oc- 
curs. More generally, strong field ionization in combina- 
tion with MFPADs may provide an informative probe of 
charge migration and charge transfer on time scales down 
to a few femtoseconds [58]. 
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APPENDIX A 

The molecular structure, population analysis, and or- 
bital energies of benzonitrile were obtained from Hartree- 
Fock calculations [47 (the molecule is oriented in the 
(F, Z) plane and the CN-end pointing in the positive Z 
axis, see Fig. |5|. The HOMO and HOMO-1 of benzoni- 
trile are shown in Fig. [71 The orbital energies, given in 
Fig. [tI agree well with the experimental ionization poten- 
tials, reported in Ref. [59], and are very close in energy 
[Ip: 9.79 eV (0.36 a.u) for HOMO; 10.06 eV (0.37 a.u) 
for HOMO-1]. Population analysis shows that while the 
HOMO-1 is localized on the aromatic ring, about 30% of 
the HOMO electron density is on the CN-end, see Fig. [7| 
The experimental 3D orientation of the molecule is such 
that the nodal planes of both HOMO and HOMO-1 or- 
bitals, coincide with the molecular plane, i.e., they have 
no electron density in the molecular plane. 

For the purpose of our Stark-shift tunneling model 
we require the permanent dipole moment and the po- 
larizabilities of benzonitrile and its cation. The nonzero 
elements of the polarizability tensor of benzonitrile are 
(in a.u.): a^^^ = 50.522, a^^^ = 87.808, a^^^ = 
125.758, ai^^ = 42.397, a^^^ = 76.282, and ai^^ = 
151.467, and the permanent dipole moment of benzoni- 
trile is 1.791 a.u. (in fair agreement with the experimen- 
tal value [60 ) and 1.822 a.u. for the cation, both pointing 
along —Z from the detector to the repeller. 
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